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1 Introduction 



The Black-Scholes equation, introduced in 1973 f2i, is a very well-known model arising in 
•^r , financial mathematics. Roughly speaking, it represents the value of an option V{S, t) which 

rj ' depends on time and on an stock price S. The BS model takes the following form: 

where r is the interest short rate and a is the volatility of the stock returns. 



If we include transaction costs in model (jl.ip then we can obtain (see [6], [7]) the following 
nonlinear version of ()1.1|) : 



2 



> 

o 

where a is an adjusted volatility. Now, if we consider the stationary version of ()1.2p then 
we obtain the following ordinary differential equation: 

S^V"{Sf +p5V"(5) + q{SV'{S) - V{S)) = 0, (1.3) 

■ ■ ~2 

^v , where p — and q — — - are constants where the volatility are included. 

jj] , 2ocr^ ba^ 

- - - In U] ! the authors studied the existence and localisation of solutions for equation (|1.3p 

in an interval [c, d], c,d > 0, with Dirichlet conditions, V(c) = 14, V{d) = Vd- To do this, 
they deal with the following equation related to ()1.3p : 



V"{x)+H{x,V{x)y{x))^Q, (1.4) 

where 

p{x)x'^ — JpixY'x^ + A.x^ q{x)\xz — y\ 
H{x,y,z) = — . 

More concretely, they construct a pair of lower and upper solutions for the differential 
equation (|1.4p with Dirichlet conditions in the case Vc <Vd and then they show that solutions 
of this equation correspond with solutions of (|1.3[) . In the present paper, we generalize those 
results in the following ways: first, we let the volatility a to depend on the stock price x, 
and so we replace constants p, q by two functions p{x),q{x)] second, we replace the Dirichlet 
conditions by functional boundary conditions, which allow us to consider a very large class 
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of problems involving equation (|1.3p . So, we deal with the following generalization of Black- 
Scholes equation with transaction costs: 

r x^V"(x)^ + p(x)x^V"(x) + q(x)(xV'(x) - V(x)) = 0, 

<^ (1.5) 

where p, q are nonnegative bounded functions in [c, d] and i?i : R x C([c, d]) — > M, i — 1,2, 
are functions which satisfy some conditions that we will state later. Under this framework, 
a large class of boundary conditions are included, namely: 

1. Dirichlet conditions: Bi{V{c), V) = V{c) - Vc, B2{V{d), V) = V{d) - Vd, 

2. Initial-integral conditions: Bi{V{c)^V) = V{c) — I k{x)V{x) dx, 

n 

3. Muhipoint conditions: B2{V{d), V) = V{d) - ^ V{xj), 

and so on. 

This paper is organized as follows: in Section 2, we consider problem (II. 5|) with Dirichlet 
conditions. So, we start from paper [^ and we generalize it by considering variable volatility 
and by dropping the condition Vc < Vd- In Section 3 we deal with problem (|1.5I) on its full 
version, that is, with functional boundary conditions. Namely, we provide a result on the 
existence of extremal solutions between lower and upper solutions, by using a generalized 
iteration with Dirichlet problems. Some examples of application are provided too. Inspired 
by [I] and [7], the trick to obtain solutions for (11.51) is to consider the following related 
problem: 



V"{x) + H{x, V{x),V'{x)) = 0, 
Bi(F(c),T/) = 0, B2{V{d),V)^0, 



(1.6) 



where 



^^, , p{x)x'^ — JpixY'x'^ + 4a;3(7(a;)|a;z — y\ 
H(x,y,z) = ;^-^ . 



(2.7) 
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2 Problem with Dirichlet conditions 

We begin by considering problem (|1.6|) with Dirichlet conditions, that is, 

(V"{x)+H{x,V{x),V'{x))^0, 
\ V{c) = Vc, V{d) = Vd. 
Notice that problem (12.71) is a particular case of (|1.6I) . with 

Bi{Vic),V) = V{c) - Vc, B2{V{d),V) = V{d) ~ Vd. 

We begin by introducing the concept of lower and upper solutions for this problem. For 
technical reasons, in Section 3 we will need these functions to have "angles" , and so we use 
the following weak definition, inspired by the idea of De Coster and Habets [1] . 
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Definition 2.1 We say that a G AC([c, d]) is a lower solution for problem \2.?^ if 

a{c) < Vc, a{d) < Vd, 
and for each xq G (c, d) one of the following conditions holds: 

1. D_a{xo) < D+a{xo); 

2. There exists an open interval Jq such that xq ^ Jq CZ {c,d), aij^ G W'^'^{Jo) and for 
almost all x d Jq we have 

a"{x) + H{x, a(x) , a (x)) > 0. 

We say that /3 G AC{[c,d]) is an upper solution for problem {2.1^ if 

m > K, I3{d) > Vd, 
and for each xq G (c, d) one of the following conditions holds: 

1. D-p{xn) > D+(3{xn); 

2. There exists an open interval Jq such that xq E Jq <Z {c,d), /3|jg G W'^'^{Jo) and for 
almost all x G Jo we have 

(3"{x)+H{x,p{x),l3'{x)) <0. 



Notice that if a, /3 are classical C^— lower and upper solutions for problem (j2.7p (as 
defined in [7] ) then they are also lower and upper solutions in the sense of Definition 12.11 

Proposition 2.1 The following holds: 



(1) If —r < — then the function 
d c 



d 



is a lower solution for problem \2. 7jj ; 
(2) Take k > such that 



(2. 



k > A / -T W max 

C"* y xe[c,d] 



^ix^-cd) + V,-Y^c 
2 d — c 



(2.9) 



where 

Then the function 
k 



max q{x). 



,. .^ 2 , f Vd-Vc k \ ^k Vd-Ve 



(2.10) 



is 



a lower solution for problem \2. 7| ); 
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(3) The function 



d— c 
is an upper solution for problem Jj^. 7p . 

Proof. 



., , Vd~Vc , dVc~cVd 
pyx) — — : X - 



(2.11) 



(1) First, notice that a'l{x) + H {x , a{x) , a' {x) — 0. On the other hand, 

Vd 
ai{c) ^ —c<Vc, ai{d) ^ d. 



(2) First, a2{c) = Vc, a2{d) = Vd- On the other hand, as (see [71 Lemma 1]) 

H{x,a{x),a'{x)) > —J — —\/\xa'{x) — a{x)\, 
we obtain 



H {x , a{x) , a' {x)) > — y^W max 

\ C'^ y xG[c,d] 



hx^^cd) + v,-y^^^c 



d ~ c 



Then, the choice of /c > imphes that 

a'^{x) + H[x,a{x),a'[x)) > 0. 
Notice that we can choose such a fc by virtue of being 



Hm —^ = +00. 



(3) It's easy to check, taking into account that 

/3"ix) + H{x, /3(x), /3'(x)) = H{x, /3(x), /3'(x)) < 
and 

/3(c) = v;, p{d)^Vd. 



D 

Theorem 2.1 If — < — then problem {2.'T^ has the extremal solutions, that is, the least 
and the greatest one, inside the interval 

[ai,/3] ={FeC2([c,d]) : ai{x) <V{x) < (i{x) for all x e [c, d]} . 



Moreover, if k > satisifies h2.S\) . then problem \2.1/^ has the extremal solutions inside the 
interval 

[a2,/3] = {F e C^{[c,d]) : a2{x) < V{x) < /3(x) for all x G [c,d]}. 
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Proof. 

Put 

Ei^{{x,y,z)e [c,d] xR2 . ai(x)<2/</3(x)}, 

and, if conditions referred above liold, 

E2^{{x,y,z)e [c,d] xR2 : a2(x) < y < /3(x)}. 
Thus defined, H is continuous in Ei and £'2. Moreover, we have that 



l-ff(a;, J/, 2:)l < \- \ -^rl^W + V Nl 

X \ X'^ X 

for all (x, y, z) G -Ei or (x, y, z) € £'2- Then, there exist 



i ^ n,ax ^ + \h-^m, B = max ^E^ (2.12) 

such that 

|if(x,2/,2)|<i + ^^ (2.13) 

for (a;,y, z) G £1 or (x,y,z) G £2. Then, function H satisfies the Nagumo condition both 
in El and £2. Using the fact that Q;i,a2 are lower solutions and /3 is an upper solution for 
problem (I2.7p . the conclusion holds by application of a well-known result (see [Sj). D 

Corollary 2.1 If k > satisifies i2.9\) . then problem \2.1^ has the extremal solutions in 
the interval [a,j3], where 

a{x) = max{ai(a;), a2(x)} 

and ai, a2, (3 are provided, respectively, by /i2.8]) . i2.10]) . i2.11]) . 

Theorem 2.2 In the conditions of Theorem \2.1[ each solution V of problem {2.1^ is convex 
and satisfies that xV'(x) — V{x) < for all x G [c,d]. 

Proof. 

The proof is the same to that done in [7] . Convexity of solutions provides from the fact 

V'ix) == -H{x, V{x), V'{x) > for ah x G [c, d]. 
On the other hand, the second assertion follows from the fact that function 

xe [c,d] I — >xV'ix)-V{x) 
is nondecreasing and 

dV'id) - V{d} < dp'{d) - V{d) ^d^-Vd = 0. 

D 

Theorem 2.3 Consider problem HI. -5]) with standard Dirichlet conditions, that is, 

{ x^Vixf +pix)x^V"{x) + q(x){xV'(x) - Vix) = 0, 

{ (2.14) 

i V{c) = K, V{d) = Vd. 
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(1) If — < — then problem {2.H^ has the extremal convex solutions in the Junctional 
interval [ai,/3], where ai and /3 are provided, respectively, by 112.8]) and V2.ll]) : 

(2) If k > satisfies i2. 9]) then problem \2.14^ has the extremal convex solutions in the 
functional interval [q:,/3], where a{x) — niax{ai(a;), 02(2;)} andai, a2, /3 are provided, 
respectively, by Hgj, JgJO)) and JitlT]) : 

(3) If —7 = — then oli is a solution of problem {2.14\ l. On the other hand, in the periodic 



cas 



e, Vc — Vd, the constant function V = Vc is a solution of \2.14-^. 



Proof. 

It suffices to check that if y is a solution of (|2.7p then V solves (|2.14[) . Indeed, by virtue of 
TheoremlU \xV'{x) - V{x)\ = ~{xV'{x) - V{x)), and so 



-pix^x"^ + \p{x)x^ - 4:q{x)x^{xV'{x) -Vix)) 
V"(x) = -H{x, V{x), V'{x)) = —^ ^--^ ^^3 ^• 

Then, V solves ((2l4l) . 

On the other hand, notice that if V is a convex solution of (|2.14p then y is a solution of 

(|2.7p too. If this is false then we have 

-p{x)x^ - \ pix)x'^ - 'iq{x)x'^{xV'{x) - Vix)) 
V"{x) = — — ^—-- ^-— — -^ < 0, 

and then V must be concave. 

To see part (3), notice that in the case — = — we have aifc) = Vc, ai(d) = Vd and 

d c 

a"(.T) = = II{x,ai{x),a[{x)). 

The periodic case is trivial too. D 

. . Vd Vc . . . 

The condition — - < — is essential and it cannot be removed at all, as we show in the 
d c 
following proposition. 

Proposition 2.2 The following conditions are equivalent: 

(1) 

d ~ c 

(2) The function ai{x) — —x is a lower solution for problem \2. 7p ; 

(3) Each solution of \2.1I^ is a convex solution of \2.14-^ . 

Proof. 

We have already shown that (1) => (2) => (3). On the other hand, if V is a solution of (|2.7p 
then V solves (I^Jil) if and only if \xV'{x) - V(x)\ = -[xV'{x) - V{x)) for all x G [c,d]. 
As the mapping x € [c,d] 1 — > xV'{x) — V{x) is nondecreasing, that condition is equivalent 
to dV'{d) - V{d) < 0, that is, V'{d) < Vd/d. As V is convex in [c,d] and ai{d) = V{d) 
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we obtain ai{x) < V(x) for all x G [c,d]. In particular, ai{c) < Vc, and so ai is a lower 
solution for problem (|2.7I) . 



Finally, if (2) holds then 

which implies (1). D 

Example 2.1 Consider problem {2.14-^ in the interval [c,d] — [2,6], with 

p{x) = 1 + a;^, q{x) — [x], 

where [•] denotes integer part, and boundary conditions V{2) = 9, 1^(6) — 1. Notice that in 
this case it is Vc > Vd- Condition h2. 9fl says that we have to get fc > such that 



k > W^Vl2fc+13, 

V 8 

so simple computations show that k ~ 12 satisfies i2.9\) . Then, function 

a{x) = 6x^ - 50x + 85 
is a lower solution for this problem. On the other hand, 

P{x) = 13-22: 
is an upper solution. 

Then, by application of Theorem \2.3\ problem {2.14-^ with p{x) = 1 + x^ , q{x) = [x], 
and boundary conditions V(2) = 9, V(6) = 1, has the extremal solutions in the functional 
interval 

max < —X, 6x^ — 50x + 85 > , 13 — 2x 
16 J 

3 Problem with functional boundary conditions 

In this section we deal with problem (I1.6P on its full expression. We begin by extending 
Definition 12. II to this case. 

Definition 3.1 We say that a £ AC{[c,d]) is a lower solution for problem U.6]) if 

Bi{a{c), a) < 0, B2{a{d), a) < 0, 

and for each xq € (c, d) one of the following conditions holds: 

1. D_a{xQ) < D+a{xQ); 

2. There exists an open interval Jq such that xq £ Jq C {c,d), aij^ € W^'^{Jo) and for 
almost all x (£ Jq we have 

a"{x) + H{x, a{x),a'{x)) > 0. 
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We say that (3 £ AC{[c,d]) is an upper solution for problem il.6]) if 
Bi(/3(c),/3)>0, B2(/3(d),/3)>0, 
and for each xq G (c, d) one of the following conditions holds: 

1. D-pixo) > D+Pixo); 

2. There exists an open interval Jq such that xq E Jo C (c, d), /S^jg G M^^'^(Jo) and for 
almost all x E Jq we have 



(3"{x) + H{x,l3{x),(3'{x)) <0. 



In the construction of our generalized iterative metliod, we will use two technical lemmas. 
First of them is the following generalization of Bolzano's theorem. 

Lemma 3.1 [H Lemma 2.3] Let a,b E R, a < b, and h : M. — > K a function satisfying 
h{a) < < h{b) and 

liminf/i(z) > h{x) > limsup/i(z) for all x G [a,b]. (3.15) 

Then there exist Ci, C2 G [a, b] such that h{ci) = = /i(c2) and if h[c) ~ for some c G [a, b] 
then ci < c < C2, that is, ci and ci are, respectively, the least zero and the greatest one of h 
in [a, b] . 

The second auxiliar result we need deals with the existence of extremal fixed points for 
nondecreasing operators defined in the space of absolutely continuous functions. 

Lemma 3.2 [8", Proposition 1.4.4] Let I C R a nonempty closed interval and [a, (3] a 
nonempty functional interval in AC{L). Assume that G : [a,/3] — > [a, /3] is a nondecreasing 
mapping and that there exists ip G L^{I, [0, +oo)) such that 

|(G7)'(a;)| < ip{x) for all 7 G [a, /3] and almost all x E I. 

Then, G has in [a, /3] the greatest, V* , and the least, 14, fixed points. Moreover, they satisfy 

K = min{F : GV < V}, V* = max{V^ : V < GV}. (3.16) 

Now we establish a new result on the existence of extremal convex solutions for problem 
(dH) 

Theorem 3.1 Assume that there exist a, jS E AG{[c,d]) which are, respectively, a lower and 
an upper solution for problem lll.6\) satisfying a{x) < /3(x) for all x E [c,d]. Put 

[a,/3] = {7 G AG{[c,d]) : a{x) < -f{x) < /3(x) for all x E [c,d]}, 

E — {y E M. : min a{x) < y < max /3{x)}, 

x(E[c,d] xG[c,d] 

and assume, moreover, that the following conditions hold: 
(Hi) For all 7 G [a, 13] and all y E E we have 

liminf i?i(z,7) > Bi{y,^) > limsupi?i(z, 7) (i = 1,2); 
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{H2) For all y G E the functions Bi{y,-) are nonincreasing in [a, 13] ii — 1,2), that is, if 
71,72 G [a,P] are such that 71(2;) < 72(2:) for all x G [c, d], then Bi{y,ji) > Bi{y,^2)- 

In these conditions, problem \1.0jl has the extremal convex solutions in [a, /3] . 

Proof. 

We define a mapping G : [a,f3] — > [a,/?] as follows: for all 7 £ [a,/?], G^ is the greatest 
convex solution in [a, (3] for the Dirichlet problem 

{ V" {x) + H{x,V{x),V' {x)) =0 for all xe \c,dV 

{ (3.17) 

i l/(c) = 7c, V{d) = 7d, 

where 7c, 7d are the greatest solutions in, respectively, [a(c),/?(c)] and [a{d),j3{d)]^ for the 
following respective algebraic equations: 

i3i(y,7) = 0, (3.18) 

B2(y,7) = 0. (3.19) 

Step 1: The mapping G is well-defined. First, by virtue of being a and j3 a lower and 
an upper solution for problem (jl.6|) and by condition {H2), we have for all 7 e [a, /3]: 

i?i(a(c),7) < B,{aic),a) < < i3i(/3(c),/3) < i?i(/3(c),7), 

S2(a(rf),7) < B2{a{d),a) < < B2(/3(d),/?) < S2(/?(d),7) 

and so condition (i?i) implies that the numbers 7c and 7^^ are well-defined, by application 
of Lemma 13.11 

On the other hand, the fact of being -fc and jd the greatest solutions of equations p.l8p - 
(|3.19P implies that a{c) < 7c, a{d) < 7^, /3(c) > 7c and /3{d) > 7^. So, a and /3 are, 
respectively, a lower and an upper solution for problem p.l7p . This guaranties that p.l7|) 
has the greatest convex solution in [a,/?]. (Notice that ()2.12|) - (j2.13p provides a Nagumo- 
type bound for H between our a and /?.) 

Step 2: G is a nondecreasing mapping. Let 71,72 G [a,/3] such that 71(2;) < 72(2;) for all 
X G [c, d] and we will show that Gji < Gj2- First, notice that 

Bi(G7i(c),72) < Si(G7i(c),7i) = 

and 

Si(/3(c),72)>Si(/3(c),/3)>0, 

so reasoning as above we obtain that G7i(c) < 672(0). In a similar way we prove that 
G7i(d) < Gj2{d). Now, assume that G71 ^ G72 and consider the function 

a{x) = max{G7i(2:), 672(2:)}. 

Thus defined, a is a lower solution for problem p.l7|) with conditions V{c) = 72c, V{d) = 72^. 
Indeed, if G71 < G72 in an interval (xi, 22) then 

a"{x) = 672(2;) =iJ(2, 072(2), (G72)'(2)) ^ H{x,a{x),a'{x)) for all 2 € (2:1,^2). 

In the same way, if G71 > G72 in an interval (2:1, ^2) then 

a"{x) — H [x , a{x) , a' (x)) for all x G (2;i,i2)- 
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On the other hand, if G7i(xo) — G^2{xq) then convexity of G^i and 6*72 imphes that 
D_d(xo) < D+a{xo). 

So, problem (|3.17l) with conditions V(c) = 72c, V{d) = 72^, has a solution in [a,/3], but this 
contradicts the fact that G72 is the greatest solution for this problem in [a,/3]. Then, we 
conclude that G^i < G72 and so G is a nondecreasing mapping. 

Step 3; G has the extremal fixed points. Let 7 e [a, /3]. For all x € [c, d] we have that 
(G-fYix) = (G7)'(c) - J^ His, G7(s), G-f'{s)) ds, 

and so 

f |/3(rf)-a(c)| |/3(c)-a(d)n , T 1 .' , 6 ^1 . 

G7 (a;) < max <^ ; , ; r + / K^ + By/s] ds, 

I d— c d — c ) Jc 

where A, B are as in (I2.12p . 

Then, by application of Lemma 13.21 we obtain that G has the extremal fixed points in [a, /3], 

say V* , K, which moreover satisfy p.l6p . 

Step 4; V* is the greatest convex solution of problem HI. 6]) in [a, /?]. First, it is clear, as 
GV* = V*, that V* is a solution of problem p.6p . Now, if V is another solution of (|1.6p 
then we have that V < GV and so (|3.16p implies that V < V* . So, V* is the greatest convex 
solution of problem (|1.6p in [a, /3]. 

To obtain the least convex solution of p.6p in [a, /3] we only have to redefine the mapping 
G in the obvious way. D 

Remark 3.1 Notice that condition (Hi) is satisfied, for example, if Bi(-,j) is continuous 
or if it has only downwards discontinuities. 

Theorem 3.2 In the conditions of Theorem \3.1l assume moreover that the function 

I ^ ■^(^) r ^1 

a\\x) = — ^7— a;, x^\c,d\, 

is a lower solution for problem lll.6\) such that a < ai. Then problem \1.5\) has the extremal 
convex solutions in [ai,/3]. If a ^ ai but (Hi) — {H2) are satisfied in [ai,l3] then problem 
lil.5\l has the extremal convex solutions in [a,/3], where a{x) = max{Q;(x), ai(a;)}. 

Proof. 

It is similar to the proof of Theorem 12.31 D 

Example 3.1 Consider problem lll.5\) in the interval [c,d] with the following boundary con- 
ditions: 

(Bi) "The initial value of the solution is one half of its mean value on the whole interval 

[c,d]"; 

{B2) "The final value of the solution has integer part 4"; 
Previous conditions can be written in the following form: 
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(Si) B,{V{c),V)^V{c)-^^J^ V{s)ds = 0; 
(B2) B2iV{d),V)^-[V{d)]+4 = 0. 

We will show that for V(d) G [4, 5) and d > 3c, a(x) — x and /3 = V(d) are a lower 

d 
solution and an upper solution for this problem. Indeed, we have —[a{d)] = —4 and 

a(c)--^— y^ a{s)ds^^^c^^j<0 
if d > 3c. On the other hand, the constant function /3 = V{d) is such that 

i3i(/3(c),/3) = iy(d)>0 

B2{/3{d),l3)^-[V{d)]+4^0. 

Finally, notice that for each V{d) G [4, 5), functions Bi and B2 satisfy conditions {Hi) — {H2) 
between a and /3. We can conclude, by application of Theorems \3.1\ and lS.Si that if d > Sc 
then problem lll.6\) with boundary conditions (Bi) — {B2) has the extremal solutions between 
a and /3 for each V{d) £ [4, 5). 
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